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Summary 


A  tracking  system  for  cancelling  turbulence -induced  beam 
wander  has  been  operated  a. id  preliminary  results  are  reported  for 
mean  irradiance  and  irradiance -fading  on  target.  The  experimental 
program  now  underway  on  wander-cancellation  and  finite -beam  effects 
is  described  in  detail. 

New  instrumental  developments  for  the  investigation  of  tur¬ 
bulence  intermittency  and  its  propagation  effects  are  described, 
including  a  multiplexed  field  array  of  microthermal  sensors  for  model- 
in,?  the  turbulence.  The  experimental  program  and  new  aspects  of 
statistical  data  interpretation  are  discussed,  with  preliminary  examples. 
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I.  Introduction 


Wo.-k  in  the  period  covered  by  this  report  inc.udes  (1)  the  can¬ 
cellation  of  atmospherically-induced  beam  wander  and  the  investigation 
of  finite-beam  effects  on  target-ir radiance  and  fading,  and  (2)  the  inves¬ 
tigation  of  the  nature  and  propagation  effects  of  the  intermittency  of 
turbulence.  Preliminary  wander -cancellation  results  are  presented, 
and  th''  ongoing  experiments  are  described  in  detail.  Illustrative  data 
related  to  intermittency  effects  are  also  described,  and  progress  on  the 
theoretical  framework,  instrumentation,  and  significance  of  this  prob¬ 
lem  are  discussed. 

II.  Beam-Wander  Cancellation  and  Finite -Beam  Effects 

The  preceding  report  on  this  program1  contains  a  detailed 
review  of  finite-beam  wander,  spread,  and  scintillation  considerations. 
The  goals  of  the  present  effort  are  (1)  to  demonstrate  the  cancellation  of 
wander,  utilizing  a  principle  of  reciprocity;2  and  (2)  to  quantitatively 
determine  finite-beam  wander  and  scintillation  fading  effects  and  relate 
these  to  theory,  much  of  which  has  only  recently  been  clarified.  1  We 
are  primarily  interested  in  the  illumination  of  a  point  on  a  target  plane, 
with  the  effects  of  tracking -out  wander,  and  the  statistics  and  spectra  of 
the  fading  mechanisms  as  a  function  of  the  relevant  parameters.  This 
includes  the  improvement  in  both  mean  irradiance  and  degree  of  fading 
obtainable  through  tracking. 

A.  Instrumentation 

3 

As  described  in  an  earlier  report,  the  laser  transmitter  con¬ 
sists  of  a  fast  servoed  tracking  system  with  a  15  cm  output  aperture.  The 
optics  incorporate  spatial  filtering  and  are  essentially  diffraction-limited. 

Quantitative  results  on  linear  and  log  fading  and  on  mean  irradi¬ 
ance  at  the  target  plane  are  obtained  with  the  system  shown  in  Figure  1. 

A  silicon  photodiode  is  dc -coupled  through  stable,  low-noise  amplifiers  to 


three  processing  circuits.  The  first  of  these  provides  an  output  propor¬ 
tional  to  the  instantaneous  logarithm  of  the  input,  while  the  second  responds 
to  the  linear  irradiance  fluctuations.  The  third  provides  a  signal  pxopor- 
ticnal  to  the  mean  irradiance.  All  three  outputs  are  recorded  on  an 

instrumentation  tape  recorder  for  statistical  and  spectral  processing  on 
a  digital  computer. 

The  target  receiver  is  normally  utilized  as  a  point  detector 
coincident  with  the  center  of  a  reflector  or  beacon  which  serves  as  a  source 
for  the  transmitter  tracking  unit.  When  this  tracking  source  is  large 
(e.g.,  several  cm),  the  small  photodetector  is  suspended  at  its  center 
(Figure  2a).  When  it  is  small,  a  45°  mirror  with  a  center  hole  is  util¬ 
ized  (Figure  2b).  These  sources  are  discussed  further  in  a  later  section. 

In  general,  beamsplitting  mirrors  are  avoided  because  scattering  off  of 
dust  and  imperfections  causes  the  detector  to  respond  to  reflector  or 
beacon  energy.  There  is  also  a  provision  for  utilizing  finite  receiver 
optics  in  order  to  s papally -smooth  the  scintillation  effects. 

The  receiver,  recording,  and  tracking  source  units  are  mounted 
in  a  mobile  van  which  may  be  used  at  various  distances  from  the  trans¬ 
mitter  (typically  on  the  order  of  1  km).  A  microthermal  probe  system  is 
used  for  the  non-optical  measurement  of  turbulence  strength,  or  equi¬ 
valently,  the  coherence  parameter  p  (see  Ref.  1). 

Data  processing  is  performed  on  a  DEC/PDP-11  computer,  which 
is  programmed  to  determine  (linear  and  log)  probability  distributions  and 
variances,  and  FFT  power  spectra. 

B.  Preliminary  Results 

The  initial  operation  of  the  tracker  electronics  is  shown  qualita¬ 
tively  in  Figure  3.  The  oscillograph  represents  the  difference  signal 
between  two  opposing  quadrants  in  the  tracking  photodetector,  and  in  the 
open-loop  case  indicates  the  apparent  angular  wandering  of  the  target. 

When  the  servo -controlled  galv.nometric  scanner  mirrors  are  activated, 
this  signal  is  substantially  eliminated;  the  remaining  high-frequency  noise 
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is  due  to  scinl  llation  and  tracker-response  as  limited  by  natural  fre¬ 
quencies  of  the  galvanometers.  The  signal  of  Figure  3a  shows  the 
characteristically  low  frequencies  (  al  Hz)  involved  in  beam  wander 
(or  image  dancing)  phenomena. 

The  effect  of  tracking  on  the  (point)  irradiance  at  the  target  is 
shown  in  Figure  4.  The  mean  irradiance  is  increased  approximately 
3  dB  and  is  stabilized  against  long-term  drift.  The  ratio  of  the  trans¬ 
mitter  diameter  (b)  to  the  coherence  radius  p  was  not  determined- 

ro  ’ 

according  to  Figure  19  of  Ref.  1,  the  improvement  in  mean  irradiance 
can  be  two  or  three  times  this  amount  when  b/p  is  properly  controlled, 
as  will  be  done  as  part  of  this  investigation. 

The  fluctuations  or  fading  shown  in  Figure  4  are  due  to  wander 
and  scintillations,  although  low-pass  filtering  has  removed  some  of  the 
latter  in  this  figure.  In  order  to  remove  scintillations  from  the  deter¬ 
mination  of  wander -fading,  the  log  amplitude  variances  with  and  without 
tracking  have  been  determined  as  a  function  of  low-pass  cutoff  frequencies. 
The  results  are  given  in  Table  I,  and  show  that,  although  scintillations 
strongly  predominated,  the  fading  related  to  wander  is  indeed  drastically 
reduced:  for  a  cutoff  frequency  of  1  Hz,  which  leaves  negligible  scintil¬ 
lation  components  but  much  of  the  wander  (Figure  3),  the  reduction  in 
variance  due  to  tracking  is  a  factor  of  nearly  20.  On  the  other  hand,  when 
scintillation  components  up  to  1  kHz  are  included,  the  reduction  in  total 
fading  is  nominal.  In  the  quantitative  experiments  in  this  investigation,  it 
will  be  essential  to  monitor  and  control  b/pQ  so  that  the  results  can  be 
related  to  theory  (Figure  20  of  Ref.  1)  and  a  substantial  reduction  of  total 
fading  demonstrated. 


Table  I.  Log  amplitude  variance  vs  cutofJ  frequency. 


Bandwidth  of  Photocurrent 

2 

<r 

no  track 

2 

(T 

track 

2 

<r 

track 

2 

(T 

no  track 

0-1  kHz 

0.  53 

0.40 

0.  75 

0-10  Hz 

0. 0047 

0. 00085 

0.  18 

0-1  Hz 

0. 0020 

0.  00011 

0.  055 

The  power  spectra  of  ir radiance  fluctuations,  with  ai.d  without 
tracking,  are  shown  in  Figure  5,  and  again  indicate  that  the  significant 
wander  components  are  below  a  few  Hz.  The  computer  program  is 
being  modified  to  expand  the  low-frequency  scale  of  such  a  plot. 

The  highly  preliminary  and  nonoptimum  results  given  above 
were  obtained  with  the  use  of  a  7  cm  retroreflector  for  a  target  beacon 
source.  As  discussed  below,  this  choice  can  lead  to  poor  results  and 
has  been  superseded.  The  preliminary  experiments  also  pointed  out 
improvements  needed  in  the  tracking  transmitter  mechanical  arrange¬ 
ment  and  electronics;  these  are  currently  being  implemented. 

C.  Experimental  Plan 

In  order  to  fully  investigate  the  potential  of  wander-tracking 
and  transmitter -aperture  control  to  improve  average  illumination  and 
fading  characteristics,  we  will  conduct  experiments  under  a  wide  range 
of  parameters,  employing  several  techniques  for  interpreting  the  data. 
Ideally,  this  will  include  all  of  the  following  elements: 

1 )  Raw  data 

The  basic  data  will  consist  of 

a)  Recording  of  log  and  linear  signals  at  the  target 

2 

b)  Microthermal  determination  of  C  or  p 

n  r  o 

c)  Recording  of  microthermal  fluctuations  A  T(t) 
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d)  Record  or  vectorscope  picture  of  (x,  y)  scanner 

input  signals  in  tracker  servo 

e)  Record  of  general  meteorological  conditions,  etc. 

2)  Variables  and  Parameters 

a)  Opcn/closed  loop  servo  (i.  e.  ,  no-track/track 

condition) 

b)  Size  of  receiver  at  target 

c)  Transmitter  size  (b) 

d)  C  2  or  p 

n  r  o 

e)  Low-pass  cutoff  frequency  of  signal  at  target 

f)  Transmitter  focus  conditions 

g)  Lateral  position  relative  to  centroid  of  target 

illumination 

h)  Range  (L) 

3)  Processing 

a)  Spectrum  of  AT,  to  monitor  presence  of  inertial 

4 

subrange  and  determine  inner  scale  (i  ) 

b)  Mean  target  irradiance  I 

2 

c)  Log  amplitude  variance  at  target  (a  ) 

d)  Power  spectrum  of  signal  at  target 

e)  Probability  distribution  of  linear  and/or  log  signal 

at  target 

4)  Interpretation 

The  major  indeperdent  parameter  in  interpreting  the 

data  is  the  ratio  b/p  (see  Ref.  i),  and  this  ratio  is 

o 

the  appropriate  abscissa  in  most  data  plots.  The 
detailed  phenomena  to  be  investigated  are  as  follows: 
a)  Effect  of  tracking  and  transmitter  size  on  mean 
irradiance. 

The  applicable  theory  for  comparison  with 
experimental  results  is  represented  by  Figs.  18, 
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19  in  Ref.  1,  and  includes  the  prediction  of 
a  substantial  increase  in  I,  using  wander- 
cancellation  and  an  appropriately-chosen 
transmitter  aperture. 

b)  Effect  of  tracking  and  transmitter  size  on  fading 

of  target  signal  due  to  wander  and  scintillation. 

The  applicable  theory  is  related  to  I'ig.  20  of 
Rif.  1.  In  addition  to  the  elimination  of  fading 
due  to  wander,  the  theory  predicts  a  substantial 
reduction  in  scintillation,  for  a  carefully  adjusted 
transmitter  with  the  proper  aperture  size.  This 
has  not  been  quantitatively  investigated  in  the 
past,  since  wander -cancellation  is  required. 

To  aid  in  these  experiments,  the  effects  cf  wander 
and  scintillation  may  be  mutually  separated  by 
reducing  the  latter  through  variable -receiver  - 
aperture  spatial -smoothing,  or  variable  low- 
pass  electronic  smoothing.  Equivalently,  the 
power  spectrum  of  fading,  including  appropriate 
l esolution  at  a  few  Hz  and  below,  may  be  utilized. 

c)  Probability  distribution  and  power  spectrum  of  wander¬ 

fading. 

The  theoretical  power  spectrum  may  be  deduced 
fr.  m  image -dancing  theory,  and  the  statistics 
are  predicted  in  the  recent  literature.  ^ 

d)  Mean-square  wander  angle  (obtained  from  item  1-d). 

The  applicable  theory  is  related  to  Fig.  21  and 
Eqs.  (24,  25,  28,  30)  of  Ref.  1. 

e)  Mean-square  instantaneous  spread  angle,  as  deduced 

from  mean  irradiance  with  wander  tracked  out 
(item  4 -a). 
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It  is  our  plan  to  cover  most  or  all  of  these  items.  In  addition, 
we  will  return  to  the  question  of  deliberate  angular -beam-dither  as  dis¬ 
cussed  in  Ref.  1.  Finally,  certain  theoretical  details  remain  to  be 
clarified,  as  die  mussed  in  that  reference  and  below. 

D.  Beacon/Reflector  Considerations 

The  subject  of  this  investigation  is  the  illumination  of  a  point  on 
a  (possibly -resolvable)  target,  and  the  development  of  sophisticated, 
image -processing  trackers  is  beyond  the  scope  of  the  effort:  the  results 
obtained  through  the  use  of  a  cooperative  reflector  or  beacon  are  valid 
and  will  be  immediately  applicable  to  such  systems. 

The  results  described  in  section  B  were  obtained  with  a  resolvable 
(7  cm)  retroreflector  at  the  target.  We  now  consider  this  to  be  a  sub¬ 
optimum  choice  of  tracking  source,  for  the  following  interrelated  reasons: 

1)  Scintillations  within  the  retro  cross-section  constitute 

a  spurious  input  to  the  tracker,  in  that  they  contribute 
to  the  apparent  centroid  motion  of  the  target. 

2)  Elements  of  beam  wander  which  are  resolvable  within 

the  retro  cross-section  will  have  an  apparent  motion 
whic~  is  two  times  that  observed  for  larger,  overall 
beam  components.  This  is  a  consequence  of  the 
lateral  displacement  of  rays  in  a  corner  cube. 

3)  The  atmospheric  bending  of  rays  may  not  be  perfectly 

correlated  over  the  entire  retro  cross-section,  so 
that  the  tracker  is  not  truly  optimizing  the  wander - 
cancellation  at  a  central  point. 

4)  The  coherent  illumination  of  a  resolvable  reflector  can 

lead  to  gross  tracker  errors  due  to  the  phenomenon  of 
"angular  scintillation.  "  The  implications  of 

reciprocity  with  relation  to  this  effect  are  not  fully 
understood,  but  it  is  clearly  desirable  to  avoid  the 


phenomenon. 


In  order  to  eliminate  these  problems,  we  can  utilize  a  sni?ll  retio. 
However,  the  resulting  SNR  in  the  simple  tracker  is  then  less  than  ideal. 
Consequently,  we  are  now  employing  a  small  beacon-laser  at  the  tarfet, 
which  does  not  invalidate  the  reciprocity  considerations  underlying  this 
work.  In  order  to  utilize  a  point  detector  for  target  irradiance  measure¬ 
ments,  the  beacon  can  be  enlarged  (Fig.  2a)  to  a  size  comparable  t''  the 
transmitter  resolution  scale  (=  1  cm)  and  defocused  sufficiently  to  elimi¬ 
nate  pointing  problems. 

E.  Scintillation  with  Large  Apertures 

As  discussed  in  Ref.  1,  the  application  of  reciprocity  theory  to 

Fried's  analysis  of  "atmospheric  Vnodulaticii  noise"  (coherent  fading)  in 

14 

an  optical  heterodyne  receiver  suggests  that  the  scintillation  of  target 
irradiance  will  grow  indefinitely  with  increasing  transmitter  aperture 
size  (b»p  ).  We  are  convinced  that  this  is  not  the  case,  which  indicates 
an  unidentified  breakdown  in  the  conditions  for  validity  of  this  particular 
analysis  at  large  b/p^.  To  support  this  contention,  we  cite  the  following 
points : 

1)  The  experiments  reported  in  Ref.  6  indicate  that  the 

scintillations  for  b  »  p^  approach  those  for  a  plane - 
wave  source  in  saturation.  17  This  will  be  further 
verified  on  the  present  program. 

2)  A  simple  consideration  of  phase -independent  oscil- 

1 8 

lators  of  identical  frequency,  representing  the 
discrete  pQ  areas  on  the  large  transmitter  (or 
reciprocal  hete rodyne -receiver)  aperture,  leads  to 
a  saturation  of  the  target  (or  heterodyne)  signal 
variance  vs  b. 

1  9 

3)  In  some  recent  Russian  theoretical  work,  which 

included  the  multiple  scattering  regime,  it  was 

shown  that  scintillations  will  saturate  (vs  b)  for 

large  values  of  the  phase  structure  function  Drt(b). 

9 
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In  addition,  the  spatial  scale  of  the  scintillation 
patches  is  predicted  to  be  the  transmitter  diffraction 
scale,  as  we  surmised  in  Ref.  6. 

F.  Cancellation  of  Beam  Wander  for  Uplink 

In  the  preceding  repoit  (Ref.  1),  it  was  stated  that  the  cancel¬ 
lation  of  beam  wander,  with  the  attendant  advantages  in  mean  illumina¬ 
tion  and  fading,  was  of  no  value  in  the  far  field  of  the  transmitter.  This 
statement  presupposes  a  constant  hjvel  ol.  turbulence  between  transmitter 
and  target  (horizontal  path),  and  must  be  modified  for  the  vertical  case. 

In  particular,  if  the  turbulence  is  confined  to  the  near  field  of  the  trans¬ 
mitter  on  an  uplink,  we  may  expect  advantages  similar  to  those  for  the 
horizontal  cafe,  with  similar  analytical  predictions  once  is  determined 
from  the  reciprocal  downlink.  This  may  be  seen  from  a  consideration  of 
this  reciprocal  downlink:  since  the  phase  structure  function  D^(b)  is  sub¬ 
stantially  equal  to  the  wave  structure  function  D(b),  the  conditions  for 
ignoring  amplitude  scintillation  effects  are  satisfied.  As  pointed  out  by 
Brown,  wander-cancellation  may  be  especially  valuable  in  an  uplink, 
due  to  the  large  turbulence  scales  which  may  be  encountered  and  which 
lead  to  geometrical  beam-bending  effects. 

III.  Turbulence  Intermittency 

The  macroscopic  intermittency  of  atmospheric  turbulence, 
which  is  observed  under  many  conditions,  gives  rise  to  large  fluctuations 
in  scintillations  and  other  propagation  effects,  and  leads  to  large  data- 
spreads  in  quantitative  experiments.  The  use  of  very-long  averaging 
times  to  eliminate  these  effects  obviously  masks  short-term  peaks  in 
the  scintillation  level,  and  hence  masks  deep  fading.  Thus,  there  is  a 
strong  motivation  to  understand  the  statistics  of  the  intermitt*  ncy  and 
its  propagation  effects;  this  has  been  further  strengthened  by  recent 
interest  in  fast,  single-shot  optical  systems,  e.  g.  ,  using  millisecond 
pulses  or  imaging,  in  which  the  atmosphere  is  effectively  frozen  so  that 
short-term  statistics  apply. 


We  define  intermittent  turbulence1  as  being  indicated  by  random 
microthermal  fluctuations  AT(t)  which  have  a  random  modulating  or 
miltiplicative -envelope  function  whose  characteristic  frequencies  are 
significantly  lower  that  the  lowe st-frequency  components  of  interest 
in  the  basic  fluct aations.  This  lower  limitation  on  "interesting"  frequen¬ 
cies  relates  to  the  outer  scale  of  the  Kolmogorov  spectrum,  or  to  the 

lowest  spatial  turbulence  scales  within  the  optical/lR  scintillation  filter 
functions. 

A  typical  microthermal  signal  showing  intermittent  turbulence 
is  given  in  Figure  6,  The  interesting  measurement  times  are  on  the 
order  of  the  envelope  fluctuation  times.  A  period  of  high  turbulence  may 
appear  for  e.g.  5  seconds  at  a  measurement  point;  if  the  wind  speed*  is 
2  m/sec,  the  turbulence  region  has  a  spatial  extent  (along  the  wind)  of 
10  m.  Although  this  region  may  be  warmer  and  hence  carry  additive 
low  frequencies  in  AT(t),  these  additive  components  do  not  affect  scintil¬ 
lation  and  are  not  of  interest;  the  multiplicative  low  frequencies  simply 
create  narrow  sidebands  around  each  of  the  higher-frequency  components 

of  AT(t),  and  do  not  noticeably  affect  the  autocorrelation  or  power  spec¬ 
trum. 

The  degree  of  meteorological  understanding  of  macroscale  tur¬ 
bulence  intermittency  is  poor,  so  that  mostly  empirical  means  are  neces¬ 
sary  for  modeling  the  turbulence  statistics.  With  the  use  of  these  models, 
and  an  understanding  of  the  relationship  to  short-term  propagation  statis¬ 
tics,  it  will  be  possible  to  generate  appropriate  data  for  inputs  to  computer 
codes  whi-h  simulate  propagation  effects. 


The  actual  rr.  .an  wind  speed  for  the  microthermal  data  of  Figs.  9 
12  was  1.  3  m/sec. 
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A.  Instrumental  Developments 

The  basic  field  site  constitutes  a  one-mile  path  over  uniform 
terrain,  with  coincident,  simultaneous  point-source  laser  beams  at 
4880  X  and  10.  6|i,  wavelengths,  and  wide -dynamic -range  receiver  elec¬ 
tronics.  In  addition,  a  microthermal  probe  system  permits  the  recording 
of  microthermal  variations  AT(t).  These  facilities  have  been  extensively 
described  in  previous  reporto. 

We  are  now  replicating  the  microthermal  probes,  in  order  to 

investigate  the  spatial-correlation  scales  of  the  short-term  turbulence 
2 

strength  (C  ).  For  later  use,  we  are  designing  a  10-probe-pair  array, 

2 

with  common  cabling  and  multiplexed  recording  of  short-term  C  from 

n 

each  pair.  The  individual  pairs  will  be  movable  and  can  be  lined  up  in 
an  "L"  configuration  parallel  and  perpendicular  to  the  wind  direction. 

The  basic  new  microthermal  system  is  shown  schematically 
in  Figure  7a.  The  probe  pair  is  differentially  ac -driven  with  a  100  kHz 
oscillator,  and  the  difference  in  probe  resistances  generates  a  net  ac 
input  to  an  amplifier.  The  output  is  synchronously  demodulated  and 
low-pass  filtered  to  remove  the  carrier,  resulting  in  a  signal  proper 
tional  to  AT(t). 

The  signal  is  then  squared,  resulting  in  an  output  propor  donal 

2 

to  AT  (t),  and  it  is  also  used  with  a  voltage -controlled  current  source  and 
capacitor  to  generate  long -term-average  signal  AT,  which  represents  the 
static  and  slowly-drifting  imbalance  of  the  probe  pair.  This  latter 
voltage,  which  may  typically  have  a  1  minute  time  constant,  is  utilized 
with  an  analog  multiplier  in  a  feedback  loop,  to  keep  the  probe  pair  in 
long-term  balance. 

Multiplexing  of  a  number  of  probe  pairs,  with  common  cabling, 
is  then  accomplished  as  shown  in  Figure  7b.  A  synchronous  multiplexing 
switch  commutates  the  probe  pair  being  demodulated,  as  well  as  the  cor- 
responding  (stored)  AT  signal  for  balance.  The  AT  (t)  signal  is  averaged 
(integrated)  over  a  single  time  window  (t  ),  and  is  then  "dumped"  to  make 
ready  for  the  next  probe  pair.  This  multiplexed  sample-mean  is  then 
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recorded  on  an  instrumentation  recorder,  along  with  the  multiplexing 
signal. 

The  demultiplexing  process  is  shown  conceptually  in  Figure  8. 
The  voltages  on  each  capacitor  represent  a  sample-hold,  short-term 
average  of  AT  (t)  (i.  e. ,  they  represent  C^ft)  )  for  the  corresponding 
probe-pair.  ^ These  signals  may  be  further  averaged  to  give  long  term 
values  of  Cn  .  or  they  may  be  cross -correlated,  as  in  the  present  appli- 
cation.  In  practice,  the  data-proces  sing  step  shown  in  Figure  8  will  be 
accomplished  on  a  digital  computer. 

If  the  number  of  probe  pairs  is  N,  the  complete  commutation 
period  Tm  is  on  the  order  of  NT  .  The  limitations  on  the  system  are  that 
11  '2TM  mu3t  be  sh°rter  than  the  variation-time -scales  of  interest  in 
Cn  (t)’  and  ,2)  *■>*  individual  time  window  T=  T^/N  must  be  somewhat 
larger  than  the  inverse  recorder  bandwidth.  Referring  to  Figure  6, 
which  is  a  typical  trace  of  intermittent  turbulence,  the  decorrelation 
times  in  Cn  (t)  are  on  the  order  of  a  few  seconds;  if  we  choose  T  =0.1 
sec  and  N  =  10.  we  have  a  recorder-bandwidth  requirement  of  onty  a 
few  hundred  Hz.  The  final  averaging,  for  noise-free  correlation  results, 
can  be  continued  for  as  long  as  the  data  record  permits;  the  effective 
averaging  period  in  real  time  T  will  be  on  the  order  to  T/N. 

u-Je  Computer  Data  Processing 


We  .nil  utilize  a  DEC/PDP-11  digital  computer  for  several  data 


processing  'unctions,  whose  applications  will  be  discussed  in  a  later  sec 
tion.  These  functions  are  currently  operational. 


1 )  Determination  of  the  probability  distributions  of 
thermal  fluctuations  (and  scintillations),  and  their 
variances.  A  number  of  successive  (1  msec)  signal- 
samples  can  be  averaged,  to  give  composite  samples; 
this  corresponds  to  variable -low-pass  filtering  or 
averaging  of  the  raw  signal,  and  can  influence  the 
probability  distribution  in  important  ways,  as  dis¬ 


cussed  below. 


2)  Direct  determination  of  higher -order  moments,  vs 

averaging  time.. 

3)  Determination  of  the  variance  of  a  sequence  of 

.  20. 

finite-observation-time  variances  (variance  noise  ), 
js  a  function  of  observation  (averaging)  time.  In  this 
case,  the  individual  samples  of  raw  data  are  not 
averaged,  but  the  observation  time  or  number  of 
samples  used  in  calculating  each  initial  variance  is 
varied,  and  then  the  "variance  of  variances"  is  cal¬ 
culated.  This  corresponds  to  squaring  before 
averaging,  and  the  technique  can  be  applied  to  higher 
moments  also. 

4)  Determination  of  the  power  spectrum  of  AT(t)  and 

AT2(t),  and  scintillations,  deluding  the  low  fre¬ 
quencies  (<1  Hz)  associated  with  the  envelope  of 
inte  r  mittenc  ie  s . 

C.  Basic  Field  Data 

The  pertinent  basic  field  data  are  as  follows. 

1)  To  model  the  intermittent  turbulent  field: 

n  I 

Record  the  multiplexed  AT  |  ,  where  the  averaging 

time  T  is  much  less  than  one  s econo,  for  probe- 
pairs  at  a  number  of  points  separated  along  lines 
parallel  and  perpendicular  to  the  wind  direction. 
Simultaneously,  measure  and  record  AT(t)  (smgle- 
probe,  no  averaging)  at  one  point. 

2  )  To  relate  the  optical/lR  scintillation  quantities  to  the 
intermittent  turbulence: 

Simultaneously  with  #1  above,  record  the  log  ampli¬ 
tude  fluctuations  I  (t)  at  the  two  wavelengths. 


In  all  cases,  conditions  with  grossly  varying  turbulence  param¬ 
eters,  such  ?.s  days  with  broken  clouds,  will  be  avoided,  since  this  does 
not  constitute  the  type  of  intermittency  of  interest.  lor  each  run,  record¬ 
ings  will  be  made  (or  interpreted)  for  as  long  as  possible  before  diurnal 
variations  are  dominant.  The  runs  should  be  made  under  a  number  of 
different  turbulence  conditions,  and  also  for  a  few  similar  conditions  in 
order  to  check  the  repeatability  of  results. 

D.  Interpretation  of  Data;  Statistical  Descriptions 

We  are  dealing  basically  with  two  random  variables:  the  log 
amplitude  optical  or  infrared  variance  f  (t),  and  the  microthermal  fluc¬ 
tuations  AT(t).  The  former  variable  is  normally  distributed,  and  the 
latter  is  approximately  (bilateral)  log  normally  distributed.21  Also,  in 
the  presence  of  intermittency,  the  latter  has  the  nature  of  a  multiplicative 
or  envelope -modulated  random  variable,  as  described  above.  In  the  non- 
saturated  scintillation  regime,  |  (ti  is  linearly  related  to  AT2  through  a 
weignted  integral  over  the  propagation  path.  Our  goals  are  to  understand 
the  statistical  nature  of  intermittency  and  its  relationship  to  averaging 
times,  scintillation  statistics,  and  experimental  data-spread. 

The  determination  of  the  statistical  properties  of  the  micro- 
thermal  fluctuations  will  hopefully  lead  to  a  meaningful  empirical  model 
for  intermittent  turbulence  near  the  ground.  In  Figure  9,  we  show  the 
probability  distribution  c/  log  |  AT(t)|  ,  for  various  averaging  times  ranging 
from  virtually  zero  averaging  (1  msec)  to  256  msec.  The  log  normal 
character  of  the  variable  is  evident  within  certain  probability  limits, 
especially  for  64  msec  averaging.  In  those  determinations,  the  total 
observation  time  was  10  minutes. 

As  the  averaging  time  increases,  corresponding  to  a  decreased 
low-pass -filter  cutoff  frequency,  we  expect  the  intermittency  to  affect 
the  statistics,  because  the  envelope  of  AT  will  show  up  as  a  slowly -varying 
mean  or  additive  low  frequency  in|  aT|.  In  fact,  for  a  sufficient  averaging 
time,  we  expect  the  statistics  of  \Zt\  -  CJt)  to  become  gaussian.  In 
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Figure  10,  we  show  the  theoretical  fourth  moment  (calculated  from 

22 

lower  moments),  for  a  log  normal  distribution,  divided  by  the  actual 
fourth  moment,  as  a  function  of  averaging  time.  The  crossover  or 
agreement  with  log  normality  occurs  at  approximate^  0.  1  sec  averaging. 
In  figure  11,  we  show  the  same  ratio  for  a  normal  distribution,  and  note 
the  approach  to  unity  approximately  20  seconds  averaging.  As  sug¬ 
gested  by  Figure  6,  this  constitutes  sufficient  averaging  of  the  inter- 
mittency  scale.  The  total  observation  time  in  these  determinations 
was  again  10  minutes. 

In  Figure  11,  we  note  an  irregularity  in  the  curve  for  averaging 

times  of  a  few  seconds.  This  corresponds  to  the  characteristic  time 

scale  of  the  intermittency.  Another  statistical  description  is  the  spread 

in  variances,  as  a  function  of  the  observation  time  for  each  variance. 

The  "variance  of  variances",  normalized  by  the  mean  variance  in  each 

case,  is  shown  vs  observation  time  in  Figure  12.  We  again  observe  an 

irregularity  or  noisiness  for  time  scales  of  a  few  seconds. 

The  examples  given  above,  although  very  preliminary,  show  tl»e 

types  of  empirical  statistics  which  can  be  used  to  model  the  intermittent 

turbulence,  and  we  discuss  these  aspects  in  more  detail  below. 

23 

Collins  has  suggested  that  intermittent  turbulence  cannot  be 

rigorously  described  by  a  Fourier  (Stieljes)  spectrum.  The  difficulty 

may  be  that  the  establishment  of  low-frequency  sidebands  on  each  basic 

component  of  tha  microth  rmal  spectrum  violates  the  condition  of  "orthog- 

24 

onal  increments".  This  may  in  turn  explain  the  formal  breakdown  of 

certain  theoretical  treatments,  such  as  those  for  averaging  time  effects. 

We  now  review  the  specific  statistical  problems  which  will  be 

addressed  in  this  investigation: 

1)  Averaging  tine  considerations 

In  the  absence  of  intermittency,  treatments  of  averaging - 

time  effects  (such  as  the  convergence  of  variance  samples 

25 

to  their  true  mean)  are  influenced  by  the  non-gaussian 
nature  of  AT(t),  although  they  will  apply  to  measurements 
of  t  (t)  or  log  amplitude  scintillation  variance.  In  the 
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1 


presence  of  intermittency,  the  treatment 
apparently  break  down  altogether,  due  to  the 
multiplicative  nature  of  the  phenomenon.  The 
resultant  micro  J  ermal  fluctuations  exhibit 
apparent  (i.  e.  ,  short-term)  nonstationarity, 
but  the  intermittency  does  not  show  up  in  the 
observed  power  spectrum  or  autocorrelation  function 
of  AT(t),  although  it  will  show  in|AT(t)|or  AT2(t). 
Formally,  the  F ourier  representation  may  be  inap¬ 
plicable,  as  discussed  above,  and  averaging  time 
considerations  may  be  largely  confined  to  empirical 
descriptions. 

2)  Measurement-sampling  statistics  (linear  regression, 

correlation,  confidence  intervals) 

We  are  interested  in  investigating  experimental 

data -spread,  including  the  iniluence  of  intermittency. 

In  the  absence  of  intermittency,  and  for  moderate 

averaging  times,  we  expect  C  J"  samples  to  be  normally 

distributed,  as  will  the  tr2  samples.  The  application 

of  standard  statistical  methods26  to  relate  empirical 

values  of  and  <r  is  straig.  liforward.  In  the 

presence  of  intermittency,  and  for  interesting  aver- 
2 

a£ing  times,  samples  will  not  be  normally  dis¬ 
tributed  (Figure  11),  and  the  same  may  be  true  of  a-2. 
Hence,  the  standard  methods  - -such  as  Lie  chi  squared 
test--require  reformulation. 

3)  Application  of  theoretical  relationships  between  C  2  and 

2  n 

(T 

In  the  nonsaturation  regime,  which  normally  applies 
at  1°  M-  wavelengths  if  not  in  the  visible,  the  two  measure 


ments  are  linearly  related  by  the  physics  of  the 
problem: 


2 

<r 


C  2(z)  f(z)  dz 
n 


(i) 


where  z  is  the  propagation-path  variable  and 

L  the  pathlength.  In  the  absence  of  intermit- 
2  2 

tency,  C  and  <r  are  normally  distributed, 
n 

and  the  above  relationship  should  be  usable  to 

2 

predict  the  statistics  of  <r  samples  (such  as 

their  spread  or  variance  vs  averaging  time) 

2 

from  statistics  of  C  samples.  It  is  clear 

n  2 

that  the  spatial  correlation  of  Cn  i^  involved, 
and  this  motivates  development  of  the  micro¬ 
probe  array  discussed  above.  In  the  presence 

of  intermittency,  the  samples  are  not  normally 

2 

distributed,  and  an  analytical  model  for 

2 

would  be  required  to  predict  complete  <r  statistics. 

This  may  be  intractable,  but  we  can  nevertheless 

,  .  ..  27 

recognize  the  relationship 

<f*2> />  -  <(Ov  c>z/.  <*> 


where  the  quantity  on  the  RHS  will  be  investigated 

with  the  microprobe  array.  Higher  moments  can 

2 

be  similarly  treated.  Hopefully,  <r  statistics  can 
be  sufficiently  well  determined  to  predict  the  per¬ 
centage  and  mean  length  of  time  that  the  short¬ 
term  variance  is  above  any  given  level,  which 
relates  to  "deep  fades, "  and  to  fast  systems  (frozen 
atmosphere). 


I 

Since  the  necessary  data  will  be  on  hand,  we  will  ai^o  determine 
2  8 

the  spectral  covariance  or  cross -wavelength-correlation  of  scintillations. 
E.  Other  Aspects 

Although  the  present  investigation  is  limited  to  scintillation,  the 
intermittency  phenomenon  will  have  important  effects  on  other  propagation 
parameters,  such  as  coherence  and  image  resolution.  An  investigation  of 
these  '  'Dies  will  be  deferred  for  the  present. 

Finally,  it  is  apparent  that  characterization  of  the  intermittency 
is  also  desirable  at  higher  altitudes. 
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VL  Figures 

1.  Simplified  schematic  diagram  of  receiver  for  measurement  of  irra- 

diance  characteristics  at  target  plane. 

2.  Mechanical  arrangements  for  tracking  source  and  detector  of  Fig.  1. 

a)  Large  tracking  scarce 

b)  Small  tracking  source 

3.  Signal  corresponding  to  angular  image  position  in  tracker.  Horizontal 

scale:  2  sec/division. 

a)  Without  tracking  servo 

b)  With  tracking  servo 

4.  Irradiance  at  target  with  and  without  tracking  <>f  atmospheric  beam- 

wander.  High  frequency  scintillations  have  been  filtered  out. 

5.  Power  spectrum  of  irradiance  fluctuations  at  taiget,  with  and  without 

tracking. 

6.  Microthermal  fluctuations  AT(t)  showing  intermittent  envelope. 

Time  scale:  5  sec/division. 

7.  Microthermal  probe  array  system. 

a)  Basic  system 

b)  Multiplexing  system 

8.  Conceptual  diagram  of  demultiplexing  process  corresponding  to  Fig.  7. 

9.  Cumulative  probability  function  for  log  |AT|  . 

a)  Averaging  times  of  1,  16,  and  256  milliseconds. 

b)  Averaging  time  of  64  milliseconds. 
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10.  Ratio  of  theoretical  fourth  moment  for  log  normal  distribution  to 

measured  fourth  moment  of  (aT|  ,  as  a  function  of  averaging 
time. 

11.  Same  as  figure  10,  for  normal  probability  distribution. 

12.  Variance  of  a  set  of  empirical  variances  of  AT,  as  a  function  of 

observation  time  for  each  empirical  variance.  The  ordinate 
is  normalized  by  the  mean  variance  for  each  averaging  time. 
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Photodetector  (Figure  1) 


FIGURE  5 


Correlator  FIGURE  8 
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